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ABSTRACT 

We have constructed a model to describe the optical emission from ultra-luminous X- 
ray sources (ULXs) . We assume a binary model with a black hole accreting matter from 
a Roche lobe filling companion star. We consider the effects of radiative transport and 
radiative equilibrium in the irradiated surfaces of both the star and a thin accretion 
disk. We have developed this model as a tool with which to positively identify the 
optical counterparts of ULXs, and subsequently derive parameters such as the black 
hole mass and the luminosity class and spectral type of the counterpart. We examine 
the dependence of the optical emission on these and other variables. We extend our 
model to examine the magnitude variation at infrared wavelengths, and we find that 
observations at these wavelengths may have more diagnostic power than in the optical. 
We apply our model to existing HST observations of the candidates for the optical 
counterpart of ULX X-7 in NGC 4559. All candidates could be consistent with an 
irradiated star alone, but we find that a number of them are too faint to fit with an 
irradiated star and disk together. Were one of these the optical counterpart to X-7, it 
would display a significant temporal variation. 
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1 INTRODUCTION 

The presence of point-like, non-nuclear and extremely lu- 
minous X-ray sources in local galaxies has been recognise d 
for some time jFabbiano fc WhiteteOollSwartz et al.l2004ft . 
These ultra-luminous X-ray sources (ULX) can have lumi- 
nosities ~ 10 40 ergs s _1 , and their true nature is yet to be 
understood. This luminosity greatly exceeds the Eddington 
luminosity Lsdd of a 10 Mq black hole (BH) if the emis- 
sion is isotropic. This has led to the suggestion that the 
accreting object could be an interm ediate mass black hole 
(IMBH) with mass 50 - 10 00 M JColbert fc Mushotzkvl 
1999: lMakishima et all EoOOL If so, these objects are a link 
between the established population of stellar mass black 
holes, and the supermassive black holes in active galactic 
nuclei. The presence of cool X-ray spectral components, and 
the times cales of rapid variabil ity in some systems (e.g. NGC 
4559 X-7. ICropper et alJl2004ft is consistent with the IMBH 
scenario. However, the existence of IMBH is still under de- 
bate, as some observations appear contradictory. For exam- 
ple, some systems require high acc retion disk temperat ures 
in model fits to the X-ray data fcbisawa et al] 12003ft . It 
has also been noted that the emission is not likely to be 
isotropic, but collimated after to a greater or lesser extent 
jKing et alJl200illKordmg et alJl2002l: lFabrikal20o1l . Alter- 



natively, it has also been argued by iBeeel man! i2002ft that 
an accretion disk dominated by radiation pressure would ex- 
hibit strong density inhomogeneities on scales much smaller 
than the disk scale height, so that an inhomogeneous accre- 
tion disk could permit escaping flux to exceed LEdd by a 
factor of up to ~ 10 — 100. 

It seems increasingly likely that the ULX population is 
heterogeneous, with evidence to support beaming with stel- 
lar m ass BHs in some cases and IMBH in others iFabbjanc 
l2004ft The QPO reported by IStrohmaver fc Mushotzkv 
(2003) is incompatible with the beaming hypothesis, for ex- 
ample, as is the association of some ULXs with diffuse Ha 
nebulae, suggesting iso tropic illumination of the interstellar 
medi um by the ULX JPakull fc Mirionil 120021 iMiller et all 
2003). In this paper we assume an IMBH interpretation 
for ULX. Furthermore, we seek to describe the very bright- 
est objects in the ULX population - specifically, those that 
have luminosities of 10 40 ergs s" 1 or greate r (e.g. in M82, 
iMatsumoto et aljEoOll : iKaaret et al] l200lft . We find that 
wind accretion from a companion star is insufficient to sup- 
ply the accretion rate required for this luminosity. We there- 
fore assume the BH accretes via Roche lobe overflow from a 
necessarily massive or giant companion. 

The IMBH interpretion depends on the X-ray emission 
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from the disk in the immediate environment of the ULX, 
which is model dependent. Moreover, the physics at these 
high accretion rates near the event horizon is far from un- 
derstood. We therefore seek an alternative channel to X-rays 
by which the nature of the ULX can be explored, using op- 
tical/infrared (optical/IR) properties. We argue that these 
properties will be strongly influenced by the proximity of 
such an intense X-ray radiation field, and this can be used 
as a diagnostic. If the ULX is in a binary system, the X-ray 
emission will modify the optical/IR characteristics of the 
companion star and accretion disk. In particular, this will 
induce intensity and colour shifts compared to normal stars, 
and these will vary at orbital periods. These both identify 
the true optical counterpart and also provide relatively di- 
rect indications of the masses of the components. 

X-ray irradiation can drive evolu tion in XRB 
flPodsi adlowski ll99ll : iRuderman et allll9 8S¥l , and cause sig- 
nificant colour and magnitude changes of the optical coun- 
terpart. This has been observed in the sub-Eddington 
regime. One example is the Her X-l system, an X-ray bi- 
nary consisting of a neutron star accreting matter from 
a non-degenerate stellar companion. T he X-ray luminosity 
is a t hird of the Eddington luminosity l|Howarth fc Wilson! 
Il983l) . and the binary period is 1.7 days. The neutron star 
accretes matter via Roche l obe overflow through an accre- 
tion disk <|Vrtilek et all200ll ). and so is a good analogy to the 
ULX systems described in this paper. The star has been ob- 
served to change spectral type fro m A to B over the binary 
period. iBahcall fc Bahcalil l)l972|) observed a B magnitude 
amplitude of 1.5 mag and interpreted this variation as a re- 
sult of X-ray heating of a late A-type star. Other authors 
interpret the variation in terms o f heating of the star and 
a tilted, precessing accret ion disk jGerend fc BovntorJll976t 

We consider a binary model for ULX to investigate the 
radiation effects. We calculate the optical/IR observables 
from the companion star and accretion disk and discuss the 
implications of the results. We apply the results to the NGC 
4559 candidates to either eliminate them as possibilities or 
constrain the parameters of any binary system of which they 
may be part. We also predict the IR properties for future 
observations. 

We describe our model in Section |5J and examine the 
results of this model in Section |2] In Section |I] we compare 
the model to observations of ULX X-7 in NGC 4559. 



2 MODEL 

The accretion rate required by the brightest ULXs exceeds 
that which could be supplied by a stellar wind. We therefore 
assume the matter is transferred onto the compact object 
through Roche lobe overflow. We will use this assumption 
to constrain the geometry of the system. We assume that the 
system is in a quasi-steady state, and the irradiated surfaces 
are in thermal, radiative, and hydrostatic equilibrium. This 
requires that the irradiated layers necessarily re-emit all of 
the radiation falling on them. 



2.1 The radiative transfer formulation 

We consider the effects of radiative transport and radiative 
equilibrium in the irradiated surface under X-ray illumina- 
tion. We consider a plane-paral lel model and ado pt the ra- 
diativ e transport formulation of lMilnel (119261) and lWu et alJ 
to describe the heated stellar surface and accretion 
disk. We detail this formulation in Appendix lAl 



2.1.1 The heated stellar surface 

The heating of any point on the stellar surface is described 
by the equations in Appendix ^ We calculated the total 
heating by dividing the surface of the star into discrete cells 
and calculating the magnitude of the effect for each cell. 
The Roche surface for a given mass ratio was calculated 
over a grid of cells. Each cell has a flat surface, the size of 
which is dependent on the distance to neighbouring points. 
The angles a and 6 (required components of the heating 
equations) can be calculated through appeal to the angle 
between the normal vector to the surface at the point in 
question, and the vector incident on the point originating at 
the BH (a) or the observer (9). The irradiation temperature 
is calculated at this point, and taken to be the temperature 
over the entire surface of the cell. 

The shape of the Roche Lobe is determined solely by 
the mass ratio q, so the angles and relative positions of the 
points are independent of the scale of the system. The sys- 
tem can therefore be easily scaled appropriately to calculate 
the heating, according to the binary separation a. This in 
turn is directly inferred from the binary period when the 
mass ratio is known. If the system is semi-detached then the 
size of the star itself is also determined by the scale of the 
system. 

We calculate the luminosity of the star by summing the 
emergent radiation in the direction of the observer for each 
cell. This includes a component as a result of irradiation, as 
well as the original stellar luminosity, which includes b oth 
limb darkening and gravity darkening ijVon ZeipeJll924Tl ef- 
fects. 



2.1.2 The heated accretion disk - Dubus et al. prescription 

We also include the emission from an irradiated accretion 
disk as an additional op tical source. We hav e considered two 
models. The first follows lPubus et alJ il999t) to describe this 
disk. There, the irradiation temperature Ti rr varies as 



4 _ r _Mc 2 



4naR 2 



(1) 



where M is the accretion rate and R is the distance from 
the accreting source. 

For accreting black holes, 

2 



r\Mc 



(2) 



where the efficiency parameter 77 ~ 0.1. If we take an X-ray 
albedo of 0.9 and assume that 77 = 0.1, then the value of C 
in Equation [Tl would be 2.57 x 10~ 3 , with the geo metry 
of a thin disk lPubus et al.lll999Ude Jong et aT1ll996D . Here 
and hereafter we use Equation [5] with r\ = 0.1 to determine 
the value of M from L x . 
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Class 


log(M/M Q ) 


log(R/R Q ) 


log(L/L Q ) 


05 V 


1.6 


1.25 


5.7 


BO V 


1.25 


0.87 


4.3 


B5 V 


0.81 


0.58 


2.9 


FO I 


1.1 


1.8 


3.9 


GO I 


1.0 


2.0 


3.8 


KO I 


1.1 


2.3 


.3.9 


MO I 


1.2 


2.7 


4.5 



Table 1. The set of parameters used in sections 151 and HI to de- 
scribe stars of different spectral type and luminosity class. Values 
are taken from lAllenl ll973l) 

Using the fact that the radiation-transfer equations are 
linear, we can use the principle of superposition to calcu- 
late the disk temperature from the combination of the irra- 
diative heating and the viscous heating in the disk in the 
absence of X-ray irradiation. We calculate the temperature 
as a function of disk radius, and then sum the flux from a 
series of blackbody annuli to describe the overall disk flux. 
We take the inner disk radius to be the last stable circular 
orbit around the BH we are describing. We take the outer 
disk radius to be the 'tidal truncation radius', beyond which 
Keplerian orbits intersect. T his is weakly dependent on the 
mass ratio JPaczvriskilll977l) but is generally taken to be be- 
tween 0.6 and 0.7 of the Roche lobe radius. We consider it 
sufficient in our model to fix the outer disk radius to be 0.6 
of the Roche lobe radius. 

2.1.3 The heated accretion disk - radiative transfer 
formulation 

The alternative description of the accretio n disk directly ap - 
plies the radiative transfer formulation of lWu et alJ j200lT) 
that we have modified for for the star (see Appendix^. We 
consider a thin disk and determine a ra dial temperature pro- 
file in absence of irradiation using the Sha kura fc Sunvaevl 
lll973l) prescrip tion. We assume th e local flare angle is given 
by h(r) oc r 9/7 foubus et alll99SQ . where h is the disk scale 
height. 

We calculate the luminosity of the disk in a manner 
identical to that of the star. We divide the disk surface into 
cells and calculate the heating effect on each as determined 
by the incident flux and the angle of incidence a. We sum 
the emergent radiation from each cell in the direction of the 
observer in order to calculate the disk luminosity. 

We will refer t o both disk models t hroughout this paper. 
The prescription of lDubus et alJ (I1999T) will be referred to as 
the first model, and the radiative transfer formulation will 
be referred to as the second. 



3 MODEL RESULTS 

In this section we examine the dependence of our model 
on various parameters. Table contains masses, radii and 
luminosities for some early-type main sequence (MS) stars 
and some supergiants. We use these stellar parameters in 
our model. We keep the primary (BH) mass as an input 
parameter, and we constrain the scale of the Roche lobes by 
setting the volume radius of the secondary lobe equal to the 



radius of the star in Tabled hence fulfilling the condition 
for Roche lobe overflow. 

We take the X-ray luminosity of the BH to be a constant 
10 40 ergs s _1 emitted isotropically. The soft X-rays are easily 
absorbed at the disk surface, whereas the hard X-rays are 
less easily absorbed but are scattered. They penetrate the 
star to greater optical depths until the photons are down 
scattered to lower energies. In our model, we parametrise the 
absorption coefficients for the hard and soft X-rays by means 
of two parameters k s and kh and denote the hard and soft 
components of the X-ray flux as Sh and S s respectively. We 
also define a hardness-ratio parameter £ = Sh/S s . We choose 
values of 2.5 and 0.01 for the two parameters k s and kh 
respectively throughout and allow the band boundary of the 
hard and soft X-ray bands as parameters to be determined. 
For an input spectrum consisting of a blackbody and a power 
law component, we find the boundary of the soft and hard 
band to be 1.5keV. We set the gravity darkening parameter 
(3 to be 0.25, representing a star with a purely radiative 
outer envelope. 

In Figure Q we show the intensity variation over the 
surface of an 05V star and disk when we take the BH mass 
to be IOMq. In Figure [2] we use the same BH mass with a 
G0I star. We use the quantity B(r) as a measure of intensity 
(equation I A9L setting r to 2/3, and we show projections of 
the star and disk in the orbital plane. The stellar maps show 
both the irradiative and darkening effects. We use the Dubus 
prescription of Section f2.1.2l to describe the disk. 

We see that the combined surface intensity is signifi- 
cantly higher than would be expected for an unirradiated 
star. There is however a noticeable difference between the 
two figures. The stellar intensity of the G0I star increases in 
the direction of the LI point, reaching a peak there. On the 
other hand, in the 05V figure the darkening effects domi- 
nate at the LI point, so that the intensity at that point is 
less than the surrounding surface. Here we are using a low 
hardness ratio of £ = 0.01, and so little flux penetrates to 
an optical depth of r = 2/3. We find that if we increase 
the hardness ratio the intensity distribution becomes simi- 
lar to that of the G0I star. As the BH mass increases, the 
separation increases, the irradiating flux decreases and the 
intensity distribution over the surface of both stars tends 
towards that shown in Figure ^b). Note that we have not 
included here any shadowing of the accretion disk on the 
stellar surface, which should magnify any darkening at the 
LI point. 

In Figure[3{a) we illustrate the change in effective lumi- 
nosity of an 05V star. We show the V band absolute mag- 
nitude against the BH mass for an unirradiated star and 
three different sets of irradiated star calculations. As well 
as our model described above, we shows the results from an 
extension to our model where we have included the effect of 
a disk which is completely opaque to the radiation incident 
on it and is hence shadowed on the stellar surface. We also 
show the effects of in addition taking into account the irra- 
diation pressure on the star - we will discuss this in section 
13.11 The disk heig ht at the outer disk radius R ou t is taken 
to be 0.2R ou t ide Jong et al.lll996f) . This shielding effect re- 
sults in a reduced stellar magnitude. In this figure we have 
set the phase angle to be zero (so the star is in superior 
conjunction) and set the inclination of the system such that 
cosi = 0.5. 
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Figure 1. The variation in intensity B(t) with r 
prescription with £ = 0.01 and a BH mass of IQMq 



= 2/3 for (b) an irradiated 05V star and (a) a disk using the iDubus et all (.1999) 
We plot projections in the orbital plane with the labelled distances in units of Rq. 
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Figure 2. The variation in intensity B(t) with r = 2/3 for (b) an irradiated G0I star and (a) a disk using the lDubus et all dlS 
prescription with £ = 0.01 and a BH mass of IOMq. We plot projections in the orbital plane with the labelled distances in units of Rq. 



Figure EI a ) shows that the heating effect on the star 
decreases with increasing BH mass, which may be counter- 
intuitive. This relationship is a consequence of constraining 
the volume radius of the secondary Roche lobe to the ra- 
dius of the undistorted star. As the mass ratio decreases, 
the Roche lobe geometry requires the binary separation a 
to increase. The result is a decrease in the amount of flux 
incident on the stellar surface. 

Figure^ b) is a sample lightcurve for the 05V star. Here 
we use a BH mass of lOOM© and an inclination such that 
cos i — 0.5. This figure shows both the ellipsoidal variation 
of an unirradiated star, as well as the combination of both 
ellipsoidal and irradiative effects. We include a third line 
showing the magnitude when the irradiated accretion disk 
is included. We will discuss the effect of the disk on the 
lightcurve in Section f3. 21 



3.1 Irradiation pressure effects 

We now consider the effect of the X-ray irradiation on the ge- 
ometrical shape of the secondary. We apply the prescription 



of lPhillips fc Podsiadlowskil J2002), which involves a modifi- 
cation of the Roche potential in which the radiation pressure 
force is parameterised using the ra tio of the radiation to the 
gravit ational force. Equation 6 of lPhillips fc Pod siadlowski 
(2002) shows how we can combine the gravitational and ra- 
diation pressure forces from the BH as a 'reduced' gravita- 
tional force equal to (1 — 5)F grav , where F grav is the grav- 
itational force. 5 is the product of a parameter dependent 
on the X-ray lumin osity and binary mass ratio (ca lculated 
using equation 25 of lPhillips fc Pod siadlows kil2002|) and the 
cosine of the angle between the surface normal and the di- 
rection of the flux vector. 

After calculating the position of each cell on the Roche 
surface in our usual way, we calculate a S value for each cell, 
replace the gravitational force in the Roche potential with 
the reduced gravitational force and recalculate the position 
of each cell. This needs to be repeated for a number of iter- 
ations in order to find the solution since the surface normal 
for each cell changes for each calculation. When the surface 
has been determined we calculate the heating at each point 
in our normal way. We take the accretion disk to be opaque 
to the radiation, so the inner Lagrangian point is shadowed 
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Figure 3. (a) The effect of irradiation on an 05V star as a function of BH mass. We take cosi = 0.5, the X-ray luminosity to be 10 4t) ergs 
s _1 emitted isotropically and the star to be at superior conjunction. The hardness ratio £ is set to 0.01. We plot lines for an unirradiated 
star, an irradiated star using our model and an irradiated star when we consider the accretion disk to be shadowed on the stellar surface. 
We also plot a line for 100.Mq ^ BH mass $C IOOO-Mq , in which we include the effect of radiation pressure 

(SectionEJ- ( b ) The V band 

absolute magnitude of the 05V star with binary orbital phase. We set the BH mass to IOOJWq, £ to 0.01 and cosi = 0.5. We show the 



magnitude variation for both an unirradiated star and for a star irradiated by a source of L x 
also show the variation when an irradiated accretion disk is included. 



10 ergs s emitted isotropically. We 



and therefore the radiation pressure will not cause the star 
to become detached from this point. 

The greater the flux incident on the stellar surface, the 
greater the pressure effect. We therefore see the most signif- 
icant distortion when the binary separation is at its lowest, 
which occurs when we use a MS star and a low BH mass. We 
illustrate this effect in Figure 21 in which we use a 150Mq 
BH and an 05V star. We show the V magnitude dependence 
on increasing BH mass in Figure |5f a). 

Note that we plot values only for BH masses of 100 - 
IOOOA/q . This is because we observe that for a BH mass of 
less than lOO Mp, the flux incident on the s urfac e is extremely 
high and the IPhillips fc Podsiadlowskl (l2002t) formulation 
becomes inappropriate to describe the stellar shape, since 
the formulation does not allow for any surface motion. In re- 
ality the external irradiation will drive circulatory currents 
in the stellar surface. A full treatment will require hydrody- 
namical motions to be considered. 



3.2 Inclusion of the accretion disk 

We now investigate further the additional flux from the ac- 
cretion disk. An increased BH mass leads to a larger binary 
separation and thus to a corresponding increase in the size of 
the accretion disk, since the outer disk radius is constrained 
by the Roche lobe size through tidal effects. The net result 
is that the disk total luminosity increases with BH mass, 
and hence compensates for the decreasing stellar total lumi- 
nosity. 

We find that the luminosity class of the irradiated star 
is the most important factor in determining which compo- 
nent dominates. To illustrate this we show in Figure |5] the 
absolute magnitude dependence on BH mass for a 05V and 
a G0I star, along with th e corre sponding disk magnitudes. 
We use the iDubus e^y l^ 99|) dis k prescription for Fig- 
ure GJa,b), and the lWu et alJ l200lf> disk model for Figure 
|5|c.d). We find the magnitude of this second disk model 
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x 

Figure 4. The variation in intensity B(t) with r = 2/3 for an ir- 
radiated 05V star with a BH mass of 150 Mp, and with disk shad- 
owing . We include the formulation of IPhillips fc Podsiadlowskil 
12002) in order to show the effect of irradiation pressure on the 
star. The solid outline shows the stellar shape when we neglect 
this effect. We plot a projection perpendicular to the orbital plane 
with the labelled distances in units of Rq ■ 



is strongly dependent on the hardness parameter £. A low 
value of £ = 0.01 produces a disk almost identical to the 
Dubus disk in terms of magnitude and colour. We use a 
value of £ = 0.1 in (c) and (d) to illustrate the effect of a 
harder X-ray spectrum. In Figure|^]we further illustrate the 
effect of varying disk hardness on the disk magnitude for 
different combinations of star and BH. 

If we examine the stellar luminosity change as a func- 
tion of BH mass in Figure 03 first, we note that while the BH 
+ MS star changes by a few tenths of a magnitude over the 
BH mass range, the BH + supergiant decreases by two mag- 
nitudes over that same range. The supergiant has a much 
larger radius, and so for a low binary separation the flux 
incident on the stellar surface will be high. However, when 
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Figure 5. V band absolute magnitudes for irradiated stars and accretion disks. We set L x = 10 4u ergs s~ L , cos(i) = 0.5 and take the 
star to be at superior conjunction. We use two different stars in our model - an Q5 type MS sta r (a,c) and a GO type supergiant (b,d). 
We plot the visu al magnitude of the star and disk against BH mass using the lDubus et alj il999l) disk formulation of Section IgTSI fa.b) 
and the lWu et alj 1200 J) disk formulation of Section 12.1,31 fed). The disk magnitude in this second formulation is very dependent on the 
hardness ratio £ of the X-ray illuminating flux (see Appendix^}- We use £ = 0.1, at which the disk magnitude is greater than for the 
first prescription. 
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Figure 6. The V band absolute magnitude of an irradiated accretion disk using the Wu et alj J200lf> disk formulation for different values 
of the hardness ratio £. We plot values for BH masses of 10, 100 and IOOOMq with an 05 type MS star (a) and a GO type supergiant 
(b). We set L x = 10 40 ergs s _1 and cos(i) = 0.5. 
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the mass ratio is decreased, this larger radius leads to a cor- 
respondingly larger binary separation than we see in the MS 
systems. 

If we now examine the disk intensity dependence on BH 
mass, we find the reverse is true. If the donor is a MS star, 
the Dubus disk (a) increases in V magnitude by more than 
2.5 magnitudes over the mass range. When we use harder 
X-ray radiation (c) the result is a more luminous disk, with 
approximately the same increase in magnitude over the mass 
range. In contrast, the Dubus disk accompanying the super- 
giant (b) increases in V magnitude by less than a magni- 
tude over the mass range (the exact increase depends on 
the method we use to determine Rout)- The disk irradiated 
by the harder X-rays (d) increases by about 1.5 magnitudes. 
These can be explained by the fact that the large supergiant 
leads to a large Roche lobe for all BH masses. Hence even a 
low BH mass results in a very bright disk, and since the tem- 
perature of the disk decreases with increasing disk radius, 
the effect of making a large disk larger still has a smaller 
effect in terms of total disk luminosity. In constrast, when 
the companion star is on the main sequence, the smaller size 
of the system at low BH masses results in a small and faint 
disk. When the BH mass is increased and the disk grows, 
the effect on its magnitude is much more significant. 

We also show on Figure [£]the V magnitude dependence 
on BH mass of the disk and star combined. It is interesting 
to note that were we actually observing an 05V system, it 
would be much easier to constrain the BH mass with the 
disk component included. The same cannot be said for the 
system with the GO supergiant. The gradient of the lumi- 
nosity change with increasing BH mass is still dictated by 
the decreasing stellar luminosity, but the curve is rendered 
shallower by the disk component. 

We now investigate the effect of the accretion disk on 
the amplitude of the lightcurve. We include a line in Figure 
E{b) showing the magnitude when the irradiated accretion 
disk is included. In our thin disk approximation the con- 
tribution of the disk will be constant for any phase, so the 
shape of the lightcurve will not be affected. The exception 
to this will be when the inclination is such that the disk is 
partially or fully eclipsed by the star. The relative amplitude 
of the lightcurve will be affected, depending on the luminos- 
ity of the disk. In Figure |3{b) for example, the amplitude 
decreases from ~ 0.1 V magnitudes to ~ 0.07. For any given 
star the magnitude of the disk will increase with BH mass, 
and so the lightcurve relative amplitude will decrease. The 
lightcurve of the star alone will decrease in amplitude with 
increasing BH mass due to the decrease in X-ray flux in- 
cident on the star - the addition of the disk will reinforce 
this. We found that the amplitude of the lightcurve for a 
GOI star with L x = 10 40 ergs s~\ cos(i) = 0.5 and £ = 0.01 
drops from ~ 1.5 Mag to ~ 0.2 Mag as we increase the BH 
mass from 10 to lOOOM©. An 05V star under the same set 
of conditions produces a lightcurve with an amplitude of 
~ 0.17 Mag for a BH mass of 1OM . This drops below 0.1 
Mag as the mass is increased to lOOM© as shown in Figure 
|^Jb), and at a mass of IOOOM0 the lightcurve is dominated 
by the ellipsoidal variation and has an amplitude of ~ 0.05 
Mag. 

Figure HJ shows the disk magnitude for a hardness ratio 
over the range of £ = 10" 4 - 10 4 using the lWu et all (1200 ll) 
formulation. We show the magnitude for a combination of 



a 10, 100 and 1000M Q BH with a 05V (a) and a GOI (b) 
star. The change in disk luminosity over this hardness range 
is large, demonstrating the importance of this factor. 

To summarise this section, we find that the stellar lu- 
minosity component is at its greatest for low BH masses and 
the disk component is at its greatest for high BH masses. If 
we consider separately a MS star, a supergiant star, a disk 
in a BH/MS system and a disk in a BH/supergiant system, 
we find the biggest changes in magnitude over the BH mass 
range occur for a supergiant star or a BH/MS disk. In gen- 
eral then, while the emission will always consist of a disk and 
a star component, the stellar component will dominate for a 
MS star / low mass BH combination, and the disk compo- 
nent will dominate in the case of a supergiant / high mass 
BH. This assumes the X-ray radiation is soft - when the 
hardness of the X-rays is increased, the contribution of the 
disk component will increase for all BH masses, and in the 
supergiant systems in particular we begin to see domination 
by the disk component over the entire mass range. 



3.3 Irradiation effects at infrared wavelengths 

We now use our model to examine the change in the pa- 
rameters of a ULX system by making predictions at other 
wavelengths, extending into the infrared. 

We have examined the magnitude change for a star and 
disk for a wavelength of 0.5 — 4.0pm, encompassing the V, 
R, I, J, H , K and L wavebands. We use the Johnson filter 
convention, with the Kron/Cousins convention for the R and 
I bands. We show plots in Figure |7| for t he stars in Table 
Q using three different BH masses and the lWu et alj (l200lT) 
disk model with £ = 0.01. As noted before, this is essentially 
interchangeable with the Dubus disk. We use a WMq BH 
in Figure a), a 100M Q BH in Figure Qb) and a 1000M Q 
BH in Figure 0c) . We incorporate shadowing of the star by 
the disk into our stellar irradiation model. 

Firstly, we see that there is a very large range in mag- 
nitude between these different systems. Secondly, we notice 
that as the mass of the BH increases, it becomes progres- 
sively harder to distinguish between different star/disk com- 
binations with a V band observation alone. Thirdly, we see 
that there is a much more clear distinction when we extend 
observations to longer wavelengths. Note that there is a clear 
separation between the MS stars and the supergiants which 
becomes more apparent as BH mass is increased. This sug- 
gests that infrared observations will have more diagnostic 
power in determining the characteristics of the ULX than 
observations at optical wavelengths. 



4 APPLICATION TO ULX X-7 IN NGC 4559 

ISoria et alJ ll2005h used HST data to study the optical envi- 
ronment of ULX X-7 in NGC 4559. They found eight possi- 
ble candidates for the ULX optical counterpart, listing the 
B, V and Ic standard magnitudes for each in table 2 of that 
paper. In this section we apply our model to this system with 
the aim of further constraining the candidate population. 

We use colour-magnitude diagrams to c ompare our 
model predictions with the observations of ISoria et alJ 
(2005). We again use the parameters of Tabled This set 
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Figure 7. The absolute magnitude of various irradiated stars and 
disks at wavelengths of 0.5 — A.Ofim. We take cos i = 0.5, the X- 
ray luminosity to be 10 40 ergs s _1 and the star to be at superior 
conjunction. We show the magnitudes (from to p to bottom) for 
BH masses of 10, 100 and 1OOOM . We use the lWu et all feOOll) 
disk prescription with £ = 0.01. 



is sufficient for us to make generalisations about the spec- 
tral type and luminosity class of the d onor star. We also add 
to this set the parameters inferred in lSoria et al.l l|200rJ ) for 
'Star No.l'; the most likely candidate for the optical coun- 
te rpart. These are determine d from the evolutionary tracks 
of iLeieune fc Schaererl ()200lf) for a non- irradiated, isolated 
star with the observed B, V and I colours. The mass is found 
to be 15 — 30Mq (we take the mass to be in the middle of this 
range), the bolometric luminosity is ~ 1.4±0.2 x 10 5 Lq and 
the effective temperature is T e // = 16000 ± 5000K. By ad- 
ditionally using these parameters we can make predictions 
for a slightly evolved star, based on a current evolutionary 
model. 

In each colour-magnitude diagram we plot the colours 
and magnitudes of seven of the eight candidate stars. 
No data are available for star no. 7, which fell on a hot 
HST/WFPC2 pixel. The error on the colour measurement 
is taken from the errors on the individual magnitude mea- 
surements, with the assumption that these errors are not 
correlated between bands. We then plot a line for each set 
of stellar parameters which we used in our model. The line 
shows the effect of varying the BH mass. The mass is varied 
from 10 — IOOOMq, and we indicate on each line where the 
mass is equal to 10, 100 and lOOOM©. 

Figure |H| compares the stars of Tabled with the obser- 
vations for systems for different inclinations, since we have 
no knowledge of the orientation of the system. We take the 
inclination to be cos(i) = 0.0 in Figure[HJa) and cos(i) = 0.5 
in Figure |HJb,c). As before we use two disk models - the first 
of which describes both the prescription of Section ^. 1.2l and 
the model of Section [2.1. 31 when £ = 0.01. The second disk 
model uses the formulation of Section ^. 1.3l but with £ = 0.1, 
giving a brighter disk. We have again taken the phase angle 
of the binary to be such that the star is at superior con- 
junction. We have examined other phase angles but we find 
that this does not affect the general results described in this 
section. We also find that the effect of including radiation 
pressure for BH masses of 100 - lOOOM© is small. We in- 
corporate shadowing of the star by the disk into our stellar 
irradiation model. 

We first examine Figure Ufa) ■ This describes the case 
where cos(z) = 0.0, so we see no accretion disk component. 
Examining the lines for the Table stars, we see a clear 
distinction between the supergiants and the MS stars. The 
MS stars occupy the left hand side of the plot only, and 
show much less change in colour and magnitude over the 
BH mass range. If we now look at the observed stars, we 
can see most of the candidate stars seem to fit into one of 
the two regimes. It is possible that either star 5 or 8 could be 
an irradiated MS star. Stars 2, 3 and 4 could all be irradiated 
supergiants (albeit for a very large BH mass). Star 1 could 
be an irradiated supergiant, or an irradiate, evolved MS star 
(track 1). 

The picture changes significantly when we consider the 
case where cos(i) = 0.5. As well as the irradiated star, there 
is a disk component present. We look first at the fainter disk 
as described in Figure IBtb). As in Figure we see that the 
inclusion of the disk results in an increase in V magnitude 
with increasing BH mass. Stars 1, 5 and 8 can still be de- 
scribed by our models.. However, the addition of the disk 
flux has a large effect on the magnitude of the supergiant 
and high-mass BH systems. The large luminosity decrease 
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Figure 8. Colour-magnitude diagram comparing observations of 
candidates for the counterpart of X-7 in NGC 4559, with model 
predictions for a range of different stars. We use the stars of Ta- 
ble Q along with a st ar with parameters derived by applying the 
evolutionary model_ojjLe je_une_fc Schaererl 1200 ll) to the observa- 
tion of star no.l ISoria^!tHlrtT2005Fl . Thishas been labelled 'Track 
#1'. For each star we plot a line showing the model result for 
a BH mass of 10 — IOOOA/q . We set the X-ray luminosity to be 
10 40 ergs s — 1 and take the star to be at superior conjunction. Fig- 
ure (a) corresponds to an inclination of cos(«) = 0.0, so there is 
no disk component pr esent in the emissio n. In (b) and (c) we set 
cos(i) = 0.5, using the Dubus et al. ( 1999) disc formulation in the 
former and the lWu et alJ 12001) formulation with £ = 0.1 in the 
latter. 



with increasing BH mass is curtailed, and the observations 
of stars 2, 3 and 4 are too faint for supergiants to be candi- 
dates. We additionally examined the case where cos(i) = 1.0, 
and found the difference between this and the cos(i) = 0.5 
case to be minor. Again, stars 1, 5 and 8 fit with predictions 
of a star and disk. Stars 2, 3 and 4 are too faint to fit with 
the predictions of the model. 

We now consider the brighter disk caused by increasing 
£ by an order of magnitude, as described in Figure Hfc). 
This disk completely dominates the optical emission - the 
only role the star plays is to constrain the disk brightness 
by determining the Roche lobe geometry. We draw the same 
conclusions from these plots as we did from Figure IHIb) - 
Stars 1, 5 and 8 fit within the predictions of these star and 
disk combinations. Stars 2, 3 and 4 do not. 

We find therefore that the observed colours of stars 2, 
3 and 4 are marginally consistent with irradiated red super- 
giants when we assume a large BH mass and set cos(i) = 0.0. 
In this case there will be temporal variations which will act 
as the candidate signature. When we move to lower incli- 
nations the increased flux from the accretion disk in our 
models indicates that the optical counterpart cannot be a 
red star (redder than B — V ~ 0.6), and so we suggest that 
star 1, 5 or 8 is more likely to be the counterpart. In this 
case, the optical observations are unable to discriminate be - 
tween different BH masses. We note that lSoria et al 
suggested without taking the effects of a disk or irradiated 
companion into account that star 1 is the most likely can- 
didate for the optical counterpart of X-7 in NGC 4559. In 
addition, taking the results of Section [3.31 into account, we 
note that the irradiated system could be very bright in the 
IR. Observations at these wavelengths could be useful in 
determining the counterpart. 



5 CONCLUSIONS 

We have constructed a model describing the heating effect 
of a ULX on a Roche lobe filling companion star. We plan 
to apply this model to the problem of positively identifying 
ULX optical counterparts, and then use observations of that 
counterpart to constrain the BH mass and the nature of the 
companion. 

Our model uses a radiative transfer formulation to ac- 
count for the X-ray nature of the incident radiation and the 
distribution of the re-radiated emission. We incorporate the 
distorted Roche lobe filling geometry of the star and account 
for the limb and gravity darkening effects. We also include 
the additional luminosity from an irradiated accretion disk, 
for two different disk models. 

We have illustrated how our model can be applied in a 
number of different ways. First, we assume that the donor 
star is filling its Roche lobe, as suggested by the high X-ray 
luminosity. We have then shown that for a given spectral 
type, the effects of irradiation decrease as the mass of the 
black hole in the system (and hence binary separation) is 
increased. Measurement of the amplitude and period of the 
lightcurve of a counterpart will allow us to determine many 
of the free parameters of the system - we have shown how 
our model incorporates both the ellipsoidal and irradiative 
effects and that we can separate out these components. We 
have discussed how the accretion disk can be the dominant 
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component in some systems and have shown that measure- 
ment of the periodicity of the binary can distinguish be- 
tween these regimes. We have shown how the importance 
of the disk is linked with the luminosity class of the com- 
panion star, and predict that observations at infrared wave- 
lengths makes it easier to distinguish between different sys- 
tems without the need for temporal observations. 

Finally, we have applied our model to a set of obser- 
vations of potential candidates for the counterpart of ULX 
X-7 in NGC 4559. While the available data is not sufficient 
to draw any firm conclusions, we see that if the observations 
contain any appreciable accretion disk component then only 
three of the possible candidates fit our model. These can- 
didates imply the donor star is an early-type MS star. Su- 
pergiant companions are only possible for a high inclination 
and a BH mass of ~ 1OOOM . 
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APPENDIX A: THE RADIATIVE TRANSFER 
FORMULATION 

We consider the effects of radiative transport in the irradi- 
ated surface under X-ray illumination. We consider a plane- 
p arallel model a nd a dopt the ra d iative transport formulation 
of lMilnd (Il926l) and IWu et all feOOll) . Milne's original for- 
mulation was for incident radiation at optical wavelengths 
and co oling via emission of optical radiation. This was mod- 
ified bvlWu" et al] l|2001) to account for incident radiation at 
X-ray wavelengths. The modification is important because 
soft X-rays will be absorbed close to the surface of the star 
by neutral and weakly ionized matter via bound-free transi- 
tions. Hard X-rays will only be attenuated at great depths 
when the matter density is significantly higher. The soft 
and hard X-ray components will subsequently have higher 
or larger opacities than for the optical radiation. The formu- 
lation is linear and therefore the principle of superposition is 
applicable. This allows us to derive the total emission using 
the irradiated and non-irradiated components. 

We take the incident radiation to be parallel beams of 
soft and hard X-rays, with effective fluxes ttS s and ttSh per 
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Figure A2. The geometry of the binary system 

unit area normal to the beams, and making an angle a to 
the normal to the stellar surface. The absorption coefficients 
of the soft and hard X-rays are k s K and k^n respectively, 
where k is the absorption coefficient of the optical radiation 
(k s > 1 and kh < 1 de fines our soft /ha rd X-ray convention 
in this study, following IWu et al1l200l|) . 

The total blackbody radiation flux is a combination of 
a component B x (r) as a result of irradiative heating by the 
incident X-rays and the component of the radiation from 
the star in the absence of irradiative heating B s (t), where 
t is the optical depth. The irradiative heating component 
B x (t) was solved in the limit of a semi infinite plane by the 
method of successive approximations and was found to be 

B x (t) — a — b s exp( — k s r sec a) — bh exp( — kur sec a) (Al) 

in the second approximation jWu et al]l200ll) . where a, b s 
and bh are constants to be determined by the boundary 
conditions. For a semi-infinite slab opaque at optical wave- 
lengths, the emergent optical/IR radiation in the direction 
8 is the Laplace transform of B x (t) 
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where A s and Ah are cosa/fc s and cos a/ kh respectively, 
and fi = cos 8. 

Here a,b s and bh are obtained by solving the radiative- 
equilibrium and radiative transfer equations for the condi- 
tions b a — » when S s — » and bh — > when Sh — > 0: 
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where the functions f 3 (a) and fh(c) are given by 
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hardness parameter £ = Sh/S s , with the total X-ray flux 
S x — Ss+Sh - By expressing B x (r) in terms of this parameter 
we obtain 



B. 
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(A8) 



As the radiative-transfer equations are linear, the local tem- 
perature stratification is given by 

T(r) = {^[B x {r) + B.(r)]j = (^Mj . (A9) 

The surface temperature of a star is effectively the temper- 
ature at an optical depth of r — 2/3. Hence, we find when 
it is viewed at a given inclination angle a, the effective tem- 
perature of the surface under irradiation is 



Teff = j— B x (2/3) +T^„ irr ^ 



1/4 



(A10) 



where T un i rr is the effective temperature in the absence of 
any irradiation. 



I - .4,, - .4;, ( Ah - ~) ln(l + k h seca) (A7) 



The hardness of the X-ray source is defined in terms of a 



